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ABSTRACT 

We have detected an Ha bow shock nebula around PSR J1741— 2054, a pulsar discovered through its GeV 
7-ray pulsations. The pulsar is only ~ 1.5" behind the leading edge of the shock. Optical spectroscopy shows 
that the nebula is non-radiative, dominated by Balmer emission. The Ha images and spectra suggest that the 
pulsar wind momentum is equatorially concentrated and implies a pulsar space velocity ^150 km s _1 , directed 
15 ± 10° out of the plane of the sky. The complex Ha profile indicates that different portions of the post-shock 
flow dominate line emission as gas moves along the nebula and provide an opportunity to study the structure 
of this unusual slow non-radiative shock under a variety of conditions. CXO ACIS observations reveal an X- 
ray PWN within this nebula, with a compact ~ 2.5" equatorial structure and a trail extending several arcmin 
behind. Together these data support a close (< 0.5 kpc) distance, a spin geometry viewed edge-on and highly 
efficient 7-ray production for this unusual, energetic pulsar. 

Subject headings: Gamma rays: stars - pulsars: individual PSR J 174 1—2054 - Shock waves 



1. INTRODUCTION 

PSR J 174 1—2054, d iscovered in a blin d pulsation search of 
a Fermi point source ( Abdo et al. 2009), is a P = 413 ms, 
characteristic age r c = 391 kyr 7-ray pulsar. The pulsar 
was then detected in archival P arkes radio observ ations and 
subsequently studied at GBT ( Cami kTet al.l 12009). showing 
it to be very faint and deeply scintillating. The radio ob- 
servations also gave a remarkably low dispersion measure 
DM = 4.7pc cm~ 3 , which for a sta ndard Galactic electron 
density model (ICordes & Laziol |2002|) implies a distance of 
« 0.38 kpc. This makes this one of the closest energetic pul- 
sars known. At this distance the low flux density observed 
at pulsar discovery, 51.4 ~ 160/iJy, also makes this the least 
luminous radio pulsar known. Indeed, since the source scintil- 
lates to undetectability at many epochs, the time average radio 
luminosity is L1.4 < 0. 025<io.4 niJy kpc 2 . The low distance 
was already suspected from the relatively high 7-ray flux (for 
the modest E = 9.4 x 10 33 erg s -1 spin-down luminosity). At 
the Z)M-estimated distance, the pulsar would emit 28% of its 
spin-down power in 7-rays, if the radiation were isotropic. 

Since PSR J 174 1-2054 is an important addition to the set 
of nearby, energetic young pulsars, further investigation of its 
energy loss and environs are needed. In particular, we wish 
to constrain its motion, true age, non-photon spindown de- 
position and interaction with its environment. For example, 
it may be a significant source of e ± injection in the solar 
neighborhood. Such injection is generally visible as a pul- 
sar wind nebula (PWN). The shape of this nebula can con- 
strain the source's distance, ISM interaction, a nd proper mo- 
tion. If the spin orientation can be measured (Ng & Romani 
2008), then we can also model the 7-ray beam shape and con- 
nect the observed 7-ray flux to its true lumin osity, predicting 
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provides the 7-ray efficiency, a critical test of pulsar models. 

We accordingly searched for a PWN, discovering an Ha 
bow shock and an X-ray PWN and trail. Like the few other 
pulsar bow shock nebulae, this is a non-radiative, B aimer- 
dominated shock. The shock velocity is relatively low, and 
Doppler broadened post-shock emission traces details of its 
flow, presenting an opportunity for the study of collision- 
less shock structure under unusual conditions, as well as con- 
straining properties of the pulsar wind. We describe initial 
observations and conclusions here. Additional discussion of 
the X-ray PWN and of the spectral properties of the pulsar 
point source are presented in Sivakoff et al.l (1201 1|) . 
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Fig. 1.— WIYN 3.6m/MiniMo Ua image of PSR J1741-2054's bow 
shock. Contours show the CXO 0.7-7 keV emission, smoothed on a 14 /; 
scale. Additional high contours smoothed at the \ " scale mark the PSR/PWN 
core. 

2. OBSERVATIONS 

2.1. Ha Imaging 

Our initial detection of the PSR J1741-2054 PWN was 
made in a 600s exposure using the WIYN 3.6m tele- 
scope and the MiniMo camera with the KP W0 12 (A = 
6566 A, AXfwhm = 16 A) narrow-band Ha filter on March 
24, 2009. An additional 4 x 300s exposure plus a continuum 
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frame were obtained the next night. Unfortunately all data 
suffered from poor (1.4— 1.8") seeing. Nevertheless, these 
data revealed a clear elliptical (13" x 20") edge-brightened 
shell. As this is close to the Galactic bulge (/ = 6.4° b = 
4.9°), the continuum frame was crowded with field stars. Fig- 
ure 1 shows the Ha image, without continuum subtraction. 
We were able to use the spectroscopic flux calibration (see 
below) to estimate the total (continuum- subtracted) Ha flux 
of the nebula as Fn a « 1.5 x 10 -14 ergcm -2 s -1 . 

We subsequently (Aug 21, 2009) obtained 3 x 600 s expo- 
sure using the EFOSC2 camera on the NTT 3.6m through the 
Ha#692 filter. While the image quality was somewhat better 
(1.1"), the wider (AXfwhm = 62 A) transmission band and 
lower peak throughput gave the Ha nebula sharper edge struc- 
ture, but lower contrast (Figure©. On this image we show the 
position of the CXO point source and the locations of the 1" 
slits used for the spectroscopic study. 




Fig. 2.— NTT 3.6m/EFOSC2 Ha image of the PSR J1741-2054 bow 
shock. Contours show the CXO 0.7-7 keV emission, smoothed on a OA" 
scale. The two rectangular boxes show the \" wide long-slit positions of the 
Keck LRIS spectroscopy. 



2.2. X-ray Imaging 

We obtained a CXO I ACIS-S observation (ObsID 11251) 
of the pulsar field on May 21, 2010 (=MJD 55337.1), with 
48.8 ks of exposure. The pulsar was placed near optimum 
focus on the backside illuminated S3 chip and the CCDs 
were operated in half-frame VFAINT mode, reading out ev- 
ery 1 .64 s to reduce pileup. A bright point source and compact 
nebula are seen, locating the pulsar to 17h 41m 57.28s, —20° 
54' 11.8" (±0.3"). The crowded optical field made it difficult 
to identify field X-ray sources. However, over the S3 half-chip 
several soft (coronal emission) sources were clearly identified 
with field stars, confirming the X-ray /optical frame tie at the 
-0.3" level. 

There is diffuse, faint X-ray emission in a PWN trail ex- 
tending some 2' at position angle PA=45 ± 5° (N through E). 
This structure contains ~ 900 ±100 counts. From 2—7" from 
the pulsar the trail contains an additional 92±14 counts. Fi- 
nally, the compact core contains ~ 3460 counts in a 2" radius 
region about the pulsar, dominated by the point source. 

The emission from these structures is consistent with a sim- 
ple absorbed power law with spectral index V w 1.6 ± 0.2 and 
absorbing column Nh ~ 1.5 ± 0.3 x 10 21 cm -2 . The point 
source clearly shows an additional soft thermal component, 
and the detailed spectral measurements of these structures are 



discussed in Sivakoff e taLl (120111) . Here we concentrate on 
the PWN morphology. 

The region around the point source appears slightly ex- 
tended. After subtracting a PSF generated with MARX, 
shifted and scaled to best match the point source compo- 
nent, the excess appears dominated by a diagonal band, « 
0.75" x 2.5", with a minor axis at PA=40 ± 5°. This struc- 
ture appears to contain ~ 400 ± 50 counts or 1 1% of the total 
counts in this region. As the minor axis PA is within 5° of the 
large scale X-ray trail axis (Figure©, we associate this struc- 
ture with the equatorial torus of the PWN. The 3:1 aspect ratio 
suggests a planar structure seen edge on, and implies that the 
pulsar spin axis is close to the plane of the sky. 




Fig. 3. — PSF- subtracted PWN core, showing extended emission (excess 
in black). If interpreted as an equatorial torus, the symmetry axis lies within 
5° of the large scale trail axis (arrow) and the Hen bow shock symmetry axis. 
The approximate location of the Hen limb is shown by the curved band. X- 
ray photons (0.7 — 7 keV) were binned to 0.25" pixels and smoothed with a 
<7 = 0.5" Gaussian after subtraction of a scaled MARX PSF. 



2.3. Nebular Spectroscopy 

Pulsar Ha bow shock nebulae are relatively rare, since the 
pulsar wind must impact a partly neutral medium. Typi- 
cally, this means that the pulsar has escaped its supernova 
remnant birthsite and is in teracting with the general ISM 
(Chatterje e & Cordesl 12002). The shock spectra, however, 
are particularly interesting as the shocked gas is advected 
downstream on a time short compared to the cooling time. 
This means that these are non-radiative shocks and are dom- 
inated byJBatoe^^missionra^ than the usual forbidden 
lines (Raymond 2001; Heng l2010h . When the up-stream neu- 
tral atoms, which are unaffected by the fields in the shock, 
drift into the shocked ISM gas they can be excited, produc- 
ing line radiation at the (cold) radial velocity of the ambient 
medium gas. For typical conditions the ratio of excitation to 
ionization rates is q ex /Qi ~ 0.25 so that we obtain an ISM- 
velocity B aimer photon for every ~ 4 neutral atoms. How- 
ever, in addition, the atoms can charge-exchange with the hot 
shocked ions. This gives rise to excited neutrals in the post- 
shock flow and a broad line component at the systemic ra- 
dial velocity of the post-shock gas. Thus spatially-resolved, 
kinematically-resolved observations of the line emission can 
provide a wealth o f information on th e local ISM and bow- 
shock kinematics (lAldcroft et al.ll2002l) . 

We were able to obtain two \" long- slit spectra of the PSR 
J1741-2054 bow shock, using the Keck I/LRIS on June 11 
2010 (MJD 55358.5). A total of 900s exposure was taken 
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with the slit at PA =50° (along the PWN symmetry axis) and 
1200s total exposure with the slit at PA =140°. Both expo- 
sures covered the pulsar position. Unfortunately, since only 
one amplifier was available in the LRIS red camera, the long 
slit was truncated. However, we managed to cover the PWN 
emission; the slit locations were accurately determined from 
recorded slit images and are shown on Figure 2. 
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Fig. 4.— Keck LRIS spectrum of the apex limb of the PSR J1741-2054 
bow shock. The extreme B aimer dominance is evident, although the strongest 
forbidden lines are weakly detected from the bright patches on the limb. 



For these observations, a 5600A dichroic separated the red 
and blue beams. The blue camera employed a 6001/mm, 
4000A blaze grating and covered 3 100-5600 A at a resolu- 
tion of A A = 0.63 A/pixel. For the bow shock exposures, the 
red camera employed the 12001/mm, 7500A grating for a res- 
olution of AA = 0.4A/pixel. With the 1" slit the achieved 
resolution as measured from the night sky lines near Ha was 
1.48 A or 68kms _1 (FWHM). The spatial resolution along 
the slit was 0.24 // /pixel. For the blue channel, we were able to 
employ standard calibrations. Unfortunately the red channel 
configuration was non-standard for the observing run and a 
flux standard was not obtained in this configuration. However, 
we were able to use flux measurements at a lower resolution 
setting along with measured grating responses to recover the 
red flux calibration with a conservative estimate of 15% accu- 
racy. The wavelength calibration was refined by reference to 
the night sky lines and line velocities hav e been converted to 
the Lo cal Standard of Rest (LSR) using the Mihalas & Binney 
(119811) constants. 

Figure H] shows the full spectrum of the bright limb in front 
of the pulsar from the PA = 50° observation. The strong 
dominance of the B aimer emission is clear, although a few of 
the brightest forbidden lines are also visible. These lines were 
detectable only in the brightest portions of the shock limb, so 
we have not been able trace this component separately. Table 
1 contains measurements of the line fluxes. 

Our B aimer measurements give Ha/Hp = 3.99 ± 0.61 and 
H7/H/3 = 0.295 ± 0.035, and we can use these flux ratios 
to estimate the extinction toward the nebula. The selective 
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TABLE 1 
-2054 Bow Shock Limb Fluxes 



Species 


Wavelenth 

A 


Flux 

10 _17 erg cm 2 s _1 




Ha 


6563 


70.2±10.5 


a 


H/3 


4861 


17.6±0.5 




H 7 


4340 


5.2±0.6 




US 


4102 


1.3±0.5 




N[II] 


6583 


4.0±0.7 


a 


N[II] 


6549 


1.0±0.4 




0[I] 


6300 


1.6±0.3 




0[I] 


6363 


0.3±0.2 




0[II] 


3727 


4.8±0.5 




S[II] 


6731 


1.0±0.2 




S[II] 


6716 


1.2±0.3 





1 Error dominated by calibration systematics 



extinction is 



E(B - V) 



2.21Log(-^f ), 

a/Pm 



-5.17Log(-^) 



where the models values (e.g. a//3 m ) depend on the shock 
populatio ns and geometry, espec ially the shock optical depth 
of Ly/3. Chevalie r et all ([1980) have estimated the B aimer 
ratios for the broad component of a radiative shock (their 
Model 2). For the optically thin (Case A) situation, they find 
at/fim = 2.9 and 7 / (3 m = 0.39. Our corresponding extinc- 
tion estimates are E(B — V) = 0.31 ± 0.15 and 0.63 ± 0.26; 
the lcr overlap range is E(B - V) = 0.37-0.46. How- 
ever if the LyP optical depth is ~ 3 — 5, trapping increases 
the expected Ha flux and the predicted broad line ratios are 
a/Pm — 4.8 and j//3 m = 0.32. The Ha values then gives 
E(B — V) <0; physical, positive values are allowed only at 
the - 1.3a level. The H7 flux gives E (B- V) = 0.18 ±0.25. 
Thus we conclude that the selective extinction is < 0.45 with 
smaller values preferred for finite thickness shocks, especially 
for the Ha measureme nt. For reference, the fitting formulae 
for extinction values in lChen etaLl ([1998) give E(B — V) = 
Ay/3.1 = 0.13 for this direction and d = 0.4 kpc. Our limit 
on the extinction can be translated to an effective Hydrogen 
column iV# « 5.6 x 10 21 E(B- V)cm- 2 < 2.5 x 10 21 cm- 2 , 
in agreement with the X-ray spectral fitting of the PWN. For 
either our conservative upper limit or the lower extinction in- 
ferred for modest optical depth shocks the neutral column is 
substantially larger than the low N e w 1.5 x 10 19 cm -2 in- 
ferred from the DM, suggesting low ionization along this line 
of sight. 

With the spectral resolution and S/N available we were able 
to resolve the Ha lines, making an initial exploration of the 
velocity structure in the post- shock flow. In Figure [5] we show 
the Ha spectra along the two slit axes, with velocity shifts 
measured relative to the LSR as determined from night sky 
lines. Emission from the approaching (blue) and receding 
(red) sides of the bow shock are clearly seen in the PA=50° 
panel. Interestingly the front of the bow shock appears to be 
dominated by blue-shifted emission. We discuss a likely in- 
terpretation of this peculiar velocity structure below. 

3. BOW SHOCK MODELING 

The spatial and kinematic information provided by these 
observations give an excellent opportunity to prob e the ge- 
ometry of the PWN outflow (lAldcroft et al.1 120021) . Very 
helpful in this study are the elegant solutions to the thin 
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mom entum-conserving bow shock provided by Wi lkinl (119961 
2000). Although these analytic results strictly only describe 
the shape of a contact discontinuity for a thin, fully mixed 
cold flow, they do provide a good approx imate shape for the 
non-radiative forward sh ock in the ISM ( Cormeron & Kapen 
1 1 9981; iBuc ciantini 2002). The general velocity structure also 
follows the numerical simulations, although the cold analytic 
model has somewhat denser flows and correspondingly higher 
velocities downstream from the bow shock apex. The charac- 
teristic standoff angle of the contact discontinuity in the thin 
shock limit is 



#o = Ro/ d ■ 



m w v w 



47rp ^ 2 



psr 



1/2 



/d=2.l"E 3i /(nv 2 7 doA) 1/2 



where for a relativistic pulsar wind the momentum flux is 
rriw^w = E/c and the numerical value is scaled to the nom- 
inal spindown flux and distance of PSR J 174 1—2054 trav- 
eling at 100^7 km s _1 through a medium of number den- 
sity n. Thus we expect a small stand-off angle for the bow 
shock, as seen in Figure 1. Numerical simulations show 
that the forward ISM shock sc ale is a factor ~ 1.3 x larger 
(Bucciantini & Bandiera 2001). 

However, Figure 1 shows that the observed Ha s hock is 
very flat at the apex; it does not fit well to the standard Wilkin 
(1996) thin shock solution. While the numerical thick bow 
shock models do show a somewhat wider angle spread down- 
stream, they do not display the very small standoff and flat- 
tened apex of the PSR J1741-2054 bow shock. Effects that 
may produce such flattened structure are underlying wind 
asymmetry and projection to the Earth line-of- sight. Wilkin 
(2000) has provided useful expressions for the structure of 
anisotropic winds. In the pulsar case, we expect that the wind 
momentum flux will be axisymmetric about the pulsar spin 
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Fig. 5. — Symmetry- (PA=50°, above) and cross- (PA=140°, below) axis 
long slit spectra of the PSR J1741— 2054 bow shock. Both pass through 
the estimated pulsar position at coordinate 0". The spectra are sky- and 
continuum- subtracted and velocity shifts (kms -1 ) are with respect to the 
LSR; residuals from bright stars appear, particularly 1.5" behind the pulsar 
(PA=50°) and 7" to the right of the pulsar (PA=140°). The symmetry-axis 
spectrum shows red- and blue-shifted components from the front and back of 
the PWN shell. The cross-axis spectrum is dominated by blue-shifted emis- 
sion near the front limb of the PWN. 



axis. In general this provides momentum deposition mis- 
aligned withjhej30^y^^ is good ev- 
idence (iJohnston et alJl2005tlNg & Romanill2008b that pulsar 
spin and proper motion axis have a tendency to be aligned. 
We therefore focus on the aligned, axisymmetric wind. 

The pulsar wind is also highly relativistic, so that the mass 
of the bow shock shell is dominated by the s wept-u p mass 
of the ISM. In the nomenclature of IWilkinl (l2QQQh this is 
a = vpsr/v w — > 0, while the aligned case is A = 0. 
We further restrict to an equatorially symmetric wind (p w = 
m w v w cx Co + c 2 cos 2 #). An isotropic wind has c 2 = while 
an equatorial (p w ex sin 2 #) wind has c 2 = —3/2. For these 
conditions the analytic solution is a surface of rotation with 
pulsar shock distance R(0) for the angle to the direction of 
the pulsar-ISM velocity with a cylindrical coordinate 



= (R/R o ) 2 sin 2 = 3(1-0 cot(9)(l-c 2 /12)+3c 2 sm 2 (9 /4. 



We can also write the tangential velocity (in the bow shock 
frame) and surface density of swept-up mass in the cold shell 
in terms of two quantities 



and 

G z = 
with 



[(1 - c 2 /12)(20 - sin2<9) + sin 3 cos<9]/4 
[(1 - c 2 /12)(l - cos2<9) + c 2 /2 (2 + cos2(9)sin 2 (9]/4 



v t = vpsr^GI + (2G Z - (S 2 ) 2 ] 1 / 2 /^ 2 
the tangential surface velocity in the shell, 

a = n fim p R (O 2 ) 1/2 v P sr/ \2v t ) 

the shell surface mass density, and nm p the mean mass per 
particle of the ambient ISM. 

This bow shock structure will be axisymmetric in some an- 
gle (j) about the pulsar velocity, which will be inclined by an- 
gle i to the Earth line-of- sight. To complete the model of the 
bow shock appearance and radial velocity structure we project 
to the plane of the sky. This places the emission at projected 
angles (C, x) wmi respect to the (pulsar) wind origin with 



and 



cos£ = cos0 cosi + sin# sim cos</ 



sinx = sin# sin0/sin£. 



The emission from each point along the bow shock is propor- 
tional to the appropriate cr(0, <j)) and the projected velocity (in 
the observer frame) is given by 



Vbs — ^t(cos7 cosi + shi7 sim cos</ 



7 psr 



with 7 = + tan" 1 [R/{dR/d0)\. 

Note that this velocity solution is for the mixed (cooled, 
zero pressure) flow of the swept up gas. In practice, nu- 
merical simulations show modest depart ures from this ana- 
lytic solution for the tangential velocity (Bucciantini 2002) 
and a detailed solution for the full line structure of the neb- 
ula should consider such effects. However, a more basic dif- 
ference lies in the fact that the non-radiative Ha arises from 
neutrals drifting into, and charge exchanging with, ions in 
the post-shock flow. Thus the Ha structure depends strongly 
on the portion of this flow which is r eached by such neutral 
atoms (Bucciantini & Bandiera 2001). In particular, immedi- 
ately behind the forward (ISM) shock we expect that the gas 
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has velocity 3/4v psr cos^, where the angle between the pulsar 
velocity and the shock normal is r] = 6 — tan -1 [(dR/dO) / R] . 
Only later does the post-shock flow converge to the tangent 
to the contact discontinuity. This effect may be seen in the 
simulations of ICormeron & Kaperl ([1998). Thus immediately 
behind the ISM shock one expects a line-of- sight velocity 

v nr = — 3/4i>p Sr cos77 cosi. 

Since cost? is always positive, this velocity always has the sign 
of — cosi. 

3.1. Comparison with the Observations 

An initial question is whether the shock created by an equa- 
torial wind can reproduce the flattened shape of this pulsar 
bow shock. One way to parametrize this shape is the ratio 
of the perpendicular half angular size 0± , measured through 
the pulsar, to the angle from the pulsar to the projected limb 
of the wind shock in the forward direction 9\\. Note that 
#|| 7^ ^osini. This ratio is very large for PSR J 1741— 2054 's 
bow shock, w 3.7 ± 0.3. We have computed this ratio for 
a number of axisymmetric, aligned wind bow shock models. 
Figure[6] shows that the projected shock limb has a ratio nearly 
independent of i for an isotropic (c2 = 0) wind, but that the 
ratio can reach the observed value for a sin 2 # (c2 = —3/2) 
wind and a pulsar motion nearly in the plane of the sky. Con- 
sidering winds that are even more equatorially concentrated 
or including the flaring effect of the finite pressure in a thick 
bow shock allows somewhat smaller i to be accommodated. 

Figure [7] shows the projected shape of an equatorial rela- 
tivistic wind emanating from the pulsar position (circle) for 
an inclination i = 80°. For comparison, the line shows 
the limb of the isotropic wind bow shock; the standoff dis- 
tance for the equatorial wind is appreciably reduced. Note 
that, since formally p w ex sin 2 #, the bow shock approaches 
the star in the forward (6 = 0) direction. Since the pul- 
sar wind likely has a jet component, a physical bow shock 
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Fig. 6. — Ratio of the perpendicular to standoff angular sizes as a function 
of i, the inclination angle of the pulsar velocity. The large ratio observed for 
PSR J 174 1—2054 is shown, suggesting an equatorially dominated wind and 
a pulsar velocity nearly in the plane of the sky. 



would be smoot h at the apex. Simil ar matching to the Ha 
limb was used by Gaen sler et al.l (120021) to argue that the wind 
of the millisecond pulsar PSR J2124— 3358 is anisotropic. For 
PSR J 174 1—2054 we have additionally been able to connect 
the shock shape with anisotropy detected in the X-ray syn- 
chrotron nebula (§2.2). 

We next check if our model can explain the features of Fig- 
ure[5j Most remarkable is the dominance of blue- shifted emis- 
sion at and in front of the pulsar. For the PA=50° slit, the 
bright knot at the apex is at v r = —34 to — 44kms _1 . The 
cross-axis slit has Ha offset to v r = —29 to — 44kms -1 , 
with no equivalent red- shifted component. This is signifi- 
cantly offset from the systemic (ISM) velocity. The lack of 
any red-shifted emission is at first surprising, given that mass 
flows around the bow shock. However, if the pulsar is moving 
out of the plane of the sky (i < 90°), we see that the prompt 
post-shock broad-line emission (from neutral atom charge ex- 
change) will show negative radial velocities from both sides of 
the nebula. This is illustrated in the top two panels of Figure 
[SI which trace the velocity structure of the prompt emission. 
Here we plot a model for i = 70° to better separate the ve- 
locity components. The i = 75 — 80° models (for a 3D space 
velocity v psr w 150 ± 50 km s _1 ) that match the shock shape 
also provide the best match to the velocity shifts in the slit 
spectra. 

The immediate post-shock layers can dominate if the neu- 
trals penetrate only partly into the shocked ISM, before suffer- 
ing nearly complete ionization This is equivalent to Case C of 
iBucciantini & Bandieral (1200 11) . In fact PSR J 174 1-2054 has 
a relatively large spindown luminosity among pulsars show- 
ing ISM Ha shocks (although below that of PSR B0740- 
28), so if the density and ionization fraction of the upstream 
medium are high, ionization may indeed be strong in the 
shocked ISM. Thus, ahead of the pulsar where the shock is 
nearly normal, neutral H may only exist in the immediate 
post-shock layer and blue shifted emission would dominate 
the observed spectrum. 

The pulsar space velocity v psr is modest and as one moves 
behind the pulsar position, the shock obliquity r] increases 
sharply and the post-shock heating drops. As ionization 
drops, neutral H may reach the bulk flow along the contact 
discontinuity (Bucciantini & Bandiera 2001 Case B). Charge 
exchange will continue and thus we expect components with 
negative radial velocity (near side) and positive radial velocity 
(far side). Indeed in Figure \5\ across the body of the shell, the 
PA=50° spectrum shows two components with velocity ex- 
trema —125 km s _1 and +25 km s -1 from the front and back 
sides of the flow. Figure [8] shows the computed radial velocity 
for this 'mixed component' . The large positive and negative 
velocities expected at the pulsar position are not seen in the 
PA=140° spectrum or in the PA=50° spectrum in front of the 
pulsar, suggesting that near the apex neutral H atoms do not 
penetrate to this 'mixed' tangential portion of the flow. The 
lower right panel shows a heuristic merged model, where the 
'mixed' component grows as 6 increases downstream. This 
smoothed model bears reasonable similarity to the observed 
2D spectrum, although it fails to reproduce the 'closed off 
shape along the symmetry axis and the red- shifted 'mixed' 
emission from the body of the nebula is less prominent in 
the data. Clearly additional effects are needed to reproduce 
the full velocity structure of this bow shock and additional 
imaging and spectroscopy will be needed to fully explain its 
origin. For example, the close-out of the symmetry axis spec- 
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Fig. 7. — Projected bow shock shape for an equatorial (c2 = —3/2) wind 
with pulsar velocity inclination i = 80° . The pulsar position is marked with 
a circle. Note that the projected bow shock limb has an appreciably flatter 
apex than that of an isotropic (63 = 0) edge-on (i = 90°) shock (line). 



trum suggests an expanding bubble structure rather than an 
open-backed bow shock. Bubble structure has been seen for 
other PWN bow shoc ks (e.g PSR B2224+65 ' s 'guit ar nebula' , 
Cordes et al. 1993); Ivan Kerkwijk & Ingld (120081) have pro- 
posed that such bubble structures are a natural consequence 
of instabilities in the post-shock flow of the relativistic pulsar 
wind. The fainter and slower red wing 'mixed' emission 5- 
10" behind the pulsar might be due to ISM inhomogeneities 
or, plausibly, instabilities in the PWN trail. 

One additional difference to the data is the lack of an ob- 
vious ISM- velocity 'narrow' line component in the measured 
spectra. In the model spectra, the component is assigned 10% 
of the broad-line flux. The limited LRIS spectral resolution 
makes it difficult to discern any such component, particularly 
in the brighter knots of the bow shock emission. However, 
we estimate from the PA=50° spectrum of the nebular body a 
broad-to-narrow intensity ratio h/I n > 3. Interestingly, this 
ratio provides a measure of T e /T p , probing electron-ion equi - 
libration in the post- shock emitting region (Raymond 2001); 
relatively large ratios and fast equilibration app ear to occur in 
low velocity shocks (Ivan Adelsberg et aLll2008|) . supporting a 
modest velocity for PSR J1741— 2054. If the narrow compo- 
nent can be isolated and kinematically resolved, it provides 
much additional information on the upstream medium's den- 
sity and on pre- shock heating of the ambient gas. 

4. GEOMETRY AND DISTANCE CONSTRAINTS 

The leading edge of the bow shock appears structured, es- 
pecially in the NTT image, so it is somewhat difficult to mea- 
sure the pulsar-apex angle and cross-axis half angle. Our best 
estimates are On w 1.5" and 8± « 5.5". We would like to 
relate these to the characteristic c ontact discontinuity standoff 
angle O . In lAldcroft et al.1 (120021) . the ISM shock/contact dis- 
continuity angle ratio was estimated as ~ 1.3 at the shock 
apex. While finite pressure effects can cause this to grow 
slightly downstream, here we assume a constant ratio. Also, 
for the C2 = — 3/2, i = 80° shape, we infer 9u = 0.5#o and 
0± = 1.8500- Together these estimates let us infer a charac- 
teristic (isotropic wind) contact discontinuity standoff angle 



O * 2.3". 

This characteristic offset implies n do. 4^7 = 0.8. As noted 
above the velocity structure suggests that « 1.5 ± 0.5, so 
we loosely constrain the product n do a « 0.4. In turn this 
suggests a plausible density n ^ 0.4 cm -3 if we adopt the 
DM distance. 

For the allowed extinction range (A(Ha) = — 1.1; i.e 
Thyp = — 5), we infer an Ha number flux from the full neb- 
ula of 7V 7 = 5 — 15 x 10 _3 cm _2 s _1 . This may be related 
to the neutral density and velocity of the upstream medium 
dRaymondll200U iGhavamian et al.1 120011) . For example the 
upstream neutrals will produce a narrow component Ha num- 
ber flux 

Kr _ q ex n H i7r(0 e ffd) 2 v _ _ 1ft _ 3in .2 „ nTn -2 Q -i 
iV 75U w q-47rd 2 — 1-oxlU uhiu^vjcvcy s 

for an excitation-to-ionization ratio q ex /qi = 0.25 and effec- 
tive full nebula angular radius 6 e ff = 1O"0io. This predicts 
< 0.3 of the observed flux of the (broad) Ha line. One can 
make a similar estimate for the flux of the broad component, 
depending on the ratio of charge-transfer to ionization rates, 
qct/q%, i.e. (q c t/qi)n H iOi v 7 = 3-10 (depending on the 
extinction). With the low vj w 1.5 and n < 1 estimated 
above, we would infer a large yield of broad H a photons 
from charge transfer intera ctions, q c t/qi > 2 — 7. Figures in 
Ivan A delsberg et all (|2008|) suggest that the broad to narrow 
ratio grows rapidly at low shock velocities, especially for op- 
tically thin shocks, although velocities as small as 150 km s _1 
are not computed. Additional modeling and, especially, an ac- 
curate measurement of the narrow line component would help 
settle whether such high efficiency conversion to Ha occurs. 

Turning to the nebula distance, we see that the relatively 
low Nh inferred from the X-ray absorption column and the 
B aimer line ratios supports a close distance for this pulsar, 
with t he HI surveys (iDickey & LockmanTl 990: Kalberla ^et al.l 
120051) showing values twice that inferred for the nebula at dis- 
tances as small as 0.5 kpc (HEASARC nH tool). With im- 
proved spectral measurements of the PWN Ha we should be 
able to fit for a precise pulsar space velocity and inclination. 
Comparison with a proper motion, measurable from CXO X- 
ray imaging (~ 7y baseline) or possibly HST imaging (few 
year baseline) would then yield a direct kinematic distance to 
the neutron star. This is particularly valuable for understand- 
ing the apparently very high 7-ray efficiency. 

PA=50, Post-shock PA=140, Post-shock 



PA=50, Mixed 




Fig. 8. — Model slit spectra for an equatorial (c2 = —3/2) wind with 
pulsar velocity inclination i = 70°. The slit axes are at 10° to the pulsar 
velocity vector and the pulsar position is marked (circle) at the systemic (nar- 
row component) velocity. Above: prompt emission from the symmetry and 
cross axes. Below: The expected velocity structure of the 'mixed' gas emis- 
sion for the PA=50° slit, and a smoothed model where the mixed gas acquires 
a neutral component, appearing downstream. 



PSR J1741-2054's Bow Shock Nebula 
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Our match to the 0\\/0± ratio suggests i « 75 ± 10°. For 
our estimated v p « 150 km s _1 this is also consistent with the 
spread of velocities in our slit spectra and with the measured 
apex blue shift - -30km s" 1 « (-3/4)cos(75°)150 kms" 1 
expected if the surrounding medium radial velocity is near 
the LSR and prompt emission dominates near the apex. If 
the pulsar spin and velocity are aligned, this implies a pulsar 
viewing angle £ > 65°. In turn this may be compared with 
the viewi ng angles inferred f or the observed 7-ray pulse pro- 
file (iRomani & WattersI 1201 til) . The two pole caustic (TPC) 
model has great difficulty producing the 7-ray pulse width 
A = 0.23 and lag from the radio 5 = 0.29. In the outer 
gap (OG) picture, 7-rays for such an old pulsar are only seen 
if viewed from near the equator, £ > 60°, although the ob- 
served pulse width and 7-radio lag are then easily achieved. 
The best-fit prefers ( « 65 — 75°, consistent with the angles 
inferred from the bow shock geometry above. 



5. CONCLUSIONS 

The discovery of an Ha bow shock associated with the 7- 
ray pulsar PSR J 174 1—2054 provides an important opportu- 
nity to constrain the geometry and momentum deposition of a 
pulsar wind. Our initial images and spectra suggest that this 
pulsar wind has a strong equatorial concentration and that the 
spin axis (and space velocity) are close to the plane of the 
sky. The spectroscopy implies a low extinction (close dis- 
tance) for the nebula, consistent with pulsar dispersion mea- 
sure estimates. Spectra were obtained for two slices across 
the nebula which are best described by a pulsar velocity of 
~ 150kms _1 directed 15 ± 10° out of the plane of the sky. 
The variation in the spectrum along the nebula's velocity axis 
suggests that different layers (prompt and then mixed) of the 
post-shock flow are emitting near the nebula apex and tail re- 



spectively. Some unmodeled features are seen and numerical 
simulation including pressure effects in the shock tail may be 
needed to explain the full velocity structure. 

The apparently equatorially concentrated wind and near 
edge-on view are consistent with inferences of PWN flow 
from observations of other pulsars and with the 7-ray emis- 
sion geometry inferred from LAT pulsations. 

Pulsar Ha bow shocks often display high symmetry sug- 
gesting that careful modeling can exploit their spatial and ve- 
locity structure to probe the dynamics of collisionless shocks 
in partly ionized media at a range of inclination angles. If, as 
suggested by our initial imaging and spectroscopy, the portion 
of the post shock flow probed by the Ha emission varies along 
the bow shock surface, then the shock of PSR J1741-2054 
provides an especially important laboratory to study the struc- 
ture of collisionless non-radiative shocks. Improved mea- 
surements could precisely determine the anisotropy of the 
relativistic pulsar wind, may constrain the degree of wind- 
velocity alignment, and can help deliver a kinematic measure- 
ment of the pulsar distance. Such studies will, however, re- 
quire higher spatial and spectral resolution data at good sensi- 
tivity. This might be acquired through intermediate resolution 
integr al field studies or Fabry-Perot Ha scans (iMann et al.l 
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